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This paper presents a simple engineering method for calculating the amount of bending 
undergone by a radio ray passing through an exponential model atmosphere. For any 
initial takeoff angle and for values of the surface refractivity ranging from 200 to 450, the 
bending angle r may be determined as a function of height above the earth's surface, using 
a few graphs and a few calculations. Indications of the accuracy of the method are given 
at the end of the paper. 



1. Introduction 

The course of a radiofrequency electromagnetic 
wave traveling through the atmosphere is altered 
by variations of the atmospheric index of refraction, 
n. These variations, due to changes in vapor 
pressure, air pressure and temperature, are extremely 
complex in detail; however, mathematical models of 
refraction can be constructed which represent an 
average picture of the variations. This paper con- 
siders an "exponential" model (the CRPL Expo- 
nential Reference Atmosphere [l] 1 ), in which the 
refractivity (n—l)X10 6 decreases exponentially with 
height, causing the radio wave to be bent away from 
its initial direction. The amount of bending is 
measured by the refraction angle, r, and is impor- 
tant in such problems as the accurate determination 
by radar of the range and height of flying objects, 
the location of extra-terrestrial radio noise sources 
in radio astronomy, and the analysis of radio com- 
munication systems. 

Mathematically, r maybe expressed m the follow- 
ing integral form [2,3] 



r = — cos 



(dn/n) [ (nr/n r ) 2 — cos 2 O ] " 

n 



1/2 



(i) 



where n is the atmospheric index of refraction, n 
is the value of n at the surface of the earth and r, 
r , and O are defined in figure 1. The CRPL Expo- 
nential Reference Atmosphere is characterized by 
an index of refraction of the form 

n=l + (n Q — l)e~ c e h 

where c e is the decay constant and h is the altitude 
above the surface of the earth. For this model at- 
mosphere eq (1) is not integrable in closed form; it 
can be expanded in series, but the resulting expres- 
sion is quite complicated for hand calculations. A 
numerical integration method has been used to com- 
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pute values of r by Bean and Thayer; these are listed 
in reference 1. This method is only practical through 
the use of a large scale computer. 

It might be noted that when O is large r may be 
calculated by a formula which is quite simple and 
very accurate [4] : 



r /no--l\ 



ctnd (l — e ° eh ) (radians). 



(2) 



However, for small S no simple expression is avail- 
able for calculations; thus, an engineering method 
was developed to provide a quick and practical 
means to obtain r in this case. This method has the 
added advantage over the tables in ref [1] in that 
the N s of r is not limited to those listed. 




Figure 1. Geometry of radio-ray refraction. 
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2. Calculation of r a 

The approximation of r is denoted by r a ; the 
formula is given by 



Ta=f(n Q ) cosd qe mn ° )+pC 
where the terms are explained as follows : 



(3) 



r a =tlie bending approximation in milliradians. 
/(^o)=read from figure 2 for a given N s . If more 
accuracy is desired, this value can be 
computed by 

«=[(?) eirK 1 )]""™-. 

n Q =the index of refraction at the earth's surface. 

&=the effective earth's radius factor. 
N s =ihe surface refractivity. 

o =the initial elevation angle expressed in milli- 
radians. 
h m =the height above the surface of the earth in 
meters. 

g=read from figure 3 for a given h m and O . 

.5= read from figure 4 or 5 for a given h m and O . 

p= a correction factor for N s which is read from 
figure 6 for a given N s . 

C= a height correction factor which is obtained 
from figure 7 for a given h m . 



f(n ) 



i — i — I — I j I. ■ . 1 r— 




Figure 2. f(n ) versus N s . 



Two examples of the computation of r a are in- 
cluded to illustrate the use of the r a formula. One 
example is for an N s of 252.9 and the second is an 
example using an N s of 404.9 in which the pO correc- 
tion factor has an effect. 

Example 1. (Calculation of r a with iV 5 < 344.5) 

Given: N s = 252.9 

0o=4O milliradians 
h m =o00 meters 

Find: r a =J(n Q ) cos 6 qe Bf(n « )+pC 

/(no) = 10.06 (fig. 2) 

cos 6>o= . 99920 
g=.0385 
B=. 00077 
Bf(n ) = . 0077 

P=0 (fig , 6) 

pC=0 
e s /(«o)+^ == 1.0077 

r a =.39 milliradians 
t=.38 milliradians 2 

Example 2. (Calculation of r a with iV 5 >344.5) 

Given: iV,=404.9 

0o=2OO milliradians 
h m =S0 meters 

Find: T a =f(n ) cos B^ge Bf{n ^ +vC 
/(%) = 24.70 
cos O -.98007 
q=. 000472 
£=.00117 



(fig- 3) 

(fig. 4) 



(fig. 2) 

(fig- 3) 
(fig. 4) 




Figure 3. q versus h m . 



2 Obtained by numerical integration. 
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Figure 1. B versus h m . 
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Figure 5. B versus h n 



Bf(n ) = . 0289 
^ = .30 
C= — .125 
p C=-m75 
Bf(n ) +pO= -.0086 

r a =.0113 milliraclians 
r = .0113 milliraclians 2 



(fig. 6) 
(fig- 7) 




Figure 6. N 8 correction factor. 
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Figure 7. Height correction factor. 

3. Derivation of r fl formula 

By plotting t versus h for many different values of 
O and n , it was decided that the simplest form that 
could be assumed for r to obtain the accuracy desired 
was 



r=f(n ) cos d e A+Bf(nQ) 



(4) 



The range of N s considered lies between 200 and 450 
since values outside this range rarely occur in actual 
practice. 

Using two values of N s , 200 and 344.5, and t's 
obtained by the numerical integration procedure of 
ref [1], a least squares fit of 
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n \/(fto) cos 6 J 



was made to obtain values of A and B for a given 
h m and O . 

Values of e A , denoted by q, were computed and 
graphed for h m and given d (see fig. 3). The B 
values were also plotted versus h m for given O (see 
figs. 4 and 5). It may be noted that q and B are not 
graphed for height values less than 10 m; at the 
upper range of height, a limit for q and B is ap- 
proached and reached for a given O - 

Using the A and B values it was found that for 
an N s greater than 344.5 a correction factor was 
needed to modify qe B /K) (or g^+6/(» )) so that the 
In (r/f(n ) cos O ) and, consequently, the r error were 
within an acceptable range. The height correction 
C was obtained from the difference between In 
(f//(^o) cos e ) and qe B f^ ) for an N s of 450 and 
plotted for graphical use (see fig. 7). Through fur- 
ther calculations and comparisons the relationship 
of N s to the height correction C was determined for 
N s greater than 344.5, and less than 450. This re- 
lationship was the basis for the N s correction p which 
was plotted versus N s (see fig. 6). Inclusion of this 
additional correction element results in an expression 
q e Bf(n )+ P c f or an jy s greater than 344.5. Since p 
equals zero when N s is less than or equal to 344.5, 

T a =f(n )cosd Q qezf^+pc 

becomes the general form for the simplified calcula- 
tion of r. 



4. Accuracy of r a Method 

In checking out the simplified calculation method 
various r a were compared with values of the CRPL 
exponential reference atmosphere r for the same h m 
and O - Of the values computed r a showed the 
smallest absolute error at the lower heights and 
smaller N s . The largest errors calculated were for 
N s of 450, the maximum absolute error being 0.53 
milliradians, with a maximum relative error of 6.8 
percent. 

Below is a table of computed values which gives 
indication of the range of absolute and percent error 
for several values of N s . It will be noted that error 
values for an N s of 200 and 344.5 were not included 
in this list since these values of N s were used for the 
least squares fit and were considered to have less 
error than the N s listed. 

Only error values for r a greater than 1.0 milli- 
radian were (arbitrarily) included in table 1. For 
r a less than 1.0 milliradian the absolute error is 
quite low, but the percent error can be high since r a 
is so small. This can give a somewhat distorted 
picture since a small r a ma} 7 have an error of only 
0.0001 milliradian and still be in error by greater 
than 3 percent. 



Table 1. — Range of error for r a >1.0 mr 



N, 


Range of abso- 
lute error 


Range of percent 
error 


252.9 
313.0 
377.2 
404.9 
450.0 


0. 0006 to 0. 0739 
to .0918 

. 0013 to . 1507 
to .2539 

.0024 to .5313 


0. 01% to 2. 11% 
to 1. 91% 

. 01% to 1. 53% 
to 1. 86% 

. 22% to 6. 79% 



5. Explanation of Symbols 



B= figures 4 and 5. 
C= height correction factor; figure 7, 
c e = decay constant; see ref [1], 
J(n )= figure 2, 

-[(IX^-"]"^ 

h = altitude above the surface of the earth, 

h m = altitude above the surface of the earth in 
meters, 

effective earth's radius factor 

= ^o 

n —r c e (n —l) 

atmospheric index of refraction 

= 1 + K-1)<T^ (CRPL Exponential 
Reference Atmosphere), 

n Q =mdex of refraction at the earth's surface 

=n(h=0), 

N s = surface refractivity 

= (7l -l)X10 6 , 

p=N s correction factor; figure 6, 





k= 



n= 



9r 
r- 



= angle at center of the earth (see fig. 1) 

^figure 3, 

= radial distance from the center of the earth, 

= distance from the center to the surface of the 
earth, 



-1/2 



= bending 

(dn/n) [(nr!n r ) 2 —cos 2 d 

^bending approximation in milliradians, 
=f(n ) cos e q e Bfin » )+pC 

- local elevation angle 

_! /Wo COS O \ / ■ G ii, 1 n 

[- -)> (using Snells law) 



=cos 



O = initial elevation or takeoff angle. 
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